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We study spatially homogeneous and isotropic solutions to the equations of motion derived from 
dilaton gravity, in the presence of a special combination of higher derivative terms in the gravitational 
action. All solutions are nonsingular. For initial conditions resembling those in the pre-big-bang 
scenario, there are solutions corresponding to a spatially flat, bouncing Universe originating in a 
dilaton-dominated contracting phase and emerging as an expanding Friedmann Universe. 
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I. INTRODUCTION 

The initial singularity is one of the outstanding prob- 
lems of current cosmological models. In standard big 
bang cosmology, the existence of the initial singular- 
ity is an inftyitable consequence of the Penrose-Hawking 
theorems. EJfl While scalar field-driven inflationary mod- 
els such as chaotic inflation B resolve many of the prob- 
lems faced by conventional cosmology, initial singular- 
ities are Sjtill generic, even when stochastic effects are 
included. H It is generally hoped that string theory may 
lead to a resolution of this problem. In an attempt to 
address the potential of string theory to remove the ini- 
tial cosmological singularity, Gasperini and Veneziaiip 
initiated a program known as pre-big-bang cosmology O 
based on the low energy effective action resulting from 
string theory. At lowest order, this is the action of dila- 
ton gravity. 



S 



1 

2^ 



1 



(1) 



where 6 is the dilaton, 



= SttG = Sttto-^, with G 



being the (4 dimensional) gravitational coupling, and 
rripi the Planck mass. The field equations of pre-big- 
bang cosmology exhibit a new symmetry, scale factor 
duality, which (in the Einstein frame) maps an expand- 
ing Friedmann-Robertson- Walker (FRW) cosmology to a 
dilaton-dominated contracting inflationary phase. This 
raises the hope that it is possible to realize a nonsin- 
gular cosmology in which the Universe starts out in a 
cold dilaton-dominated contracting phase, goes through 
a bounce and then emerges as an expanding FRW Uni- 
verse.0 

Unfortunately, it has been shown that the two branches 
of pre-big-bang cosmology cannot Jpe-, smoothly con- 
nected within the tree-level action. Ert2l The contract- 
ing dilaton-dominated branch has a future singularity. 



*See Ref. for a recent review of pre-big-bang cosmolog 



whereas the expanding branch emerges from a past sin- 
gularity. One-loop effects-,in superstring cosmology can 
regulate the singularity and smoothly connect a con- 
tracting phase to an expanding phase, at least in the 
presence of spatial curvature. Refs ||l4-P7|l describe other 
attempts to regulate the singularities of pre-big-bang cos- 
mology. 

A natural approach to resolving the singularity prob- 
lem of general relativity is to consider an effective the- 
ory of gravity which contains higher order terms, in ad- 
dition to the Ricci scalar of the Einstein action. This 
approach is well motivated, since we expect that any ef- 
fective action for classical gravity obtained from string 
theory, quantum gravity, or by integrating out matter 
fields, will contain higher derivative terms. Lastly, in the 
case of string theory, t-duality provides further evidence 
that physical quantities will remain finite at all times. 
Thus, it is extremely natural to consider higher derivative 
effective gravity theories when studying the properties of 
space-time at large curvatures. 

The various possible extensions to classical general rel- 
ativity all hold out the hope of a fully nonsingular theory 
of gravity, and therefore a nonsingular cosmology. How- 
ever, this promise has yet to be realized in an effective 
theory of gravity that is rigorously derived from a well 
motivated model of Planck-scale physics. In this paper 
we approach the problem from a different perspective, 
and ask instead whether it is possible to derive a higher 
order theory of gravity which will produce a cosmology 
that captures the qualitative features of the nonsingu- 
lar cosmological evolution envisaged by the pre-big-bang 
scenario. 

An approach to explicitly constructing an effective 
gravitational action which insures that physical invari- 
ants are always finite is given in Refs [^,0 . The result- 
ing action includes a particular combination of quadratic 
invariants of the Riemann tensor added to the usual 
Einstein-Hilbert action for gravity. This term forces all 
solutions of the equations of motion to approach de Sit- 
ter space-time at high curvature, and therefore renders 
them nonsingular. The model thus obtained is a specific 
higher-derivative gravity theory. 
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We should mention that a further difficulty with a per- 
turbative approach to the cosmological singularity prob- 
lem is that the perturbation expansion is expected to 
break down at energy scales lower than the scale at which 
the singularity will be smeared out. Refs also ad- 

dress the problem of ensuring that there is maximum 
allowable curvature in a given cosmological model, us- 
ing a similar technique to the one employed to eliminate 
singularities. At this point, however, we have chosen to 
focus our attention on the singularity problem alone, and 
have not imposed a limit on the curvature, beyond that 
implied by the removal of singularities which ensures that 
for any given solution the curvature will be bounded. 

The simplest way of achieving a nonsingular cosmology 
is to add an invariant I2 to the action, with the property 
that /2 = is true if and only if the spacetime is a de 
Sitter space. By coupling I2 into the gravitational action 
via a Lagrange multiplier field V' with a potential chosen 
to ensure that /2 — > at large curvatures we can require 
that all solutions approach de Sitter space at large curva- 
ture, thereby removing the singularity. For homogeneous 
and isotropic spacetimes, a choice for I2 which satisfies 
this condition is 

h = y^4i?^.i?^--i?2 (2) 

The simplest way of writing the resulting action is 

S{g^,,,il:) = J d^x^{R + + (3) 

where V{ip) is a function chosen such that the action 
has the correct (for V' — ^ 0) Einsteinian low curvature 
limit, whereas for jV'l — > 00 the constraint equation forces 

In this paper, we investigate the consequences of 
adding the same higher derivative terms to the action 
of pre-big-bang cosmology, and examine whether it can 
eliminate the singularities and produce a smooth bounce 
connecting a contracting phase to an expanding Universe. 
Our main result is that it is possible to find a potential 
V{ili) for the Lagrange multipliers which ensures that all 
cosmological solutions of the extended action for dila- 
ton gravity are nonsingular. Moreover, we show that 
there is a class of solutions corresponding to a contract- 
ing Universe smoothly connected to an expanding FRW 
phase via a bounce. Furthermore, this happens even in 
the absence of spatial curvature. Our model, therefore, 
constitutes a successful implementation of the goals of 
pre-big-bang cosmology. 

II. ACTION AND EQUATIONS OF MOTION 

Our starting point is the action (|^) for dilaton grav- 
ity (written in the Einstein frame) to which we add the 
higher derivative term given by I2, in analogy to what 
was done in the absence of the dilaton in Refs. [ p9| , p0[ : 



(4) 

For the moment, we allow a general coupling between I2 
and the dilaton. Minimal coupling corresponds to setting 
the constant 7 equal to zero. The constant c rescales the 
Lagrange multiplier field tp, and will be chosen to simplify 
the equations of motion. 

Restricted to a homogeneous and isotropic metric of 
the form 

ds^ ^ dt^ - a(t)^( ^-—dr^+r^dn^), (5) 

where dfl^ is the metric on S"^, the equations of motion 
resulting from (|^) become 

4> + + 7«/'e'^*yi2 - ij) = 0, 

TT _ k _ e-''^ dV 
~ ^ cVT2 04: ' 

6-^ + 6i/2 _ ^ _ v{tp) = 
ce^*VT2 (3i/2^ - ^V- + + 701^)) , (6) 

where dots denote derivatives with respect to time, t. 

Since our initial goal is to construct a spatially flat, 
bouncing Universe, we set the curvature constant fc = 0. 
We will, for the moment, consider minimal coupling of 
to I2, which fixes 7 = 0. To eliminate useless constant 
coefficients in the equations of motion, it is convenient 
to choose c\/l2 = 1. The resulting equations of motion 
become 

^ = -3i7V +GH - -|(ix' + V^(V')), 
H = -F'(V), 

X - -3ffx, (7) 

with X — s-iid a prime (/) signifying derivatives with 
respect to ip. 

We next turn to a discussion of the criteria which the 
potential V{'ip) must satisfy. At small curvatures, the 
terms in the action (^ which depend on must be neg- 
ligible compared to the usual terms of dilaton gravity. 
This is ensured by demanding 

~ V' IV-I - (8) 

as the region of small {ipl will correspond to the low cur- 
vature domain. O In order to implement the limiting 
curvature hypothesis, the invariant I2 must tend to zero, 
and thus the metric g^j/ will tend to a de Sitter metric 
at large curvatures, i.e. for [ipl 00. From the varia- 
tional equation with respect to 'ip, it is obvious that this 
requires 

V{ip) const IV-I 00. (9) 
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Conditions (||) and (|) are the same as those employed 
in Refs but they do not fully constrain the po- 

tential. In order to obtain a bouncing solution in the 
absence of spatial curvature it is necessary to add a third 
criterion: the equations must allow a configuration with 
H = and V' 7^ 0. From the equation of motion for ip in 
it follows that V{^p) must become negative, assum- 
ing that it is positive for small {tpl. Let ^l>b denote the 
nontrivial zero of V{'ip): 



Vi^bb) = 0. 



(10) 



In the absence of the dilaton, ipb will correspond to the 
value of ip at the bounce. In the presence of (j), the value 
of I '01 at the bounce will depend on x and will be larger 
than \il>b\- 

A simple potential which satisfies the conditions (||), 
(i) and (Hol) is 



1 + 




(11) 



(which is 



Note that the potential used in Refs. 
slightly simpler) does not satisfy condition ([10| 



III. PHASE DIAGRAM OF SOLUTIONS IN THE 
ABSENCE OF THE DILATON 

The conditions, (||) - (pX|), on the potential, V{'4'), dis- 
cussed in the previous section are necessary but not suf- 
ficient to obtain a nonsingular cosmology. These condi- 
tions ensure that all solutions which approach large val- 
ues of Itpl are nonsingular, but the possibility of geodesic 
incompleteness for solutions which always remain within 
the small j-i/'l region remains to be studied. In this sec- 
tion we determine the phase plane (?/;, H) of our model 
in the absence of the dilaton. We use analytical and nu- 
merical methods to study the trajectories of solutions of 
(0) in the phase plane and thus explicitly demonstrate 
the absence of singularities. Our model therefore yields 
a further example of a higher derivative gravity theory 
without cosmological singularities. In addition, and un- 
like the model of Refs. |2^j3^, we shall show that our 
theory admits spatially fiat bouncing solutions. 

There are several special points and curves on the 
phase plane {ip,H). First, the point {i^,H) = (0,0) cor- 
responds to Minkowski space-time. The potential ^(V') 
vanishes at this point, but it also vanishes at the points 



i'b = ±4. 



(12) 



As discussed earlier, the phase plane points (V'bjO) cor- 
respond to bouncing points of cosmological trajectories. 

The derivative of V{ip), and hence H, vanishes for the 
values 



(13) 



The phase plane lines {ipd, H) are therefore lines along 
which H = 0. 

To demonstrate that the point H) = (4, 0) is in 
fact a bounce, we expand the ■0 equation of motion near 
H = 0, which yields 



Hi; ^ -V. 



(14) 



Thus, as we cross the H = axis, the sign of ip changes. 
Contracting solutions with 2 < -0 < 4 have tp > and 
approach the point (4, 0) in finite time since H is posi- 
tive and does not tend to zero, and emerge for iJ > as 
expanding trajectories with decreasing curvature (since 
■0 < 0). The trajectories in the phase plane are symmet- 
ric about the H = axis, except that the time arrows 
are reversed. 

Next, we expand the equations near the origin of the 
phase plane and obtain 



dH 



1 
2H 



(0 



6IP 



(15) 



from which we see immediately the existence of critical 
lines located at 



0c(i?) = ±V6H. 



(16) 



Focusing on the contracting solutions, the equation of 
motion for H reduces to 



H 



(17) 



Consequently, trajectories which lie above the critical line 
have H < and (from Eq. (p^) ip > 0. These trajec- 
tories thus are directed towards the line = 2 where H 
changes sign. Provided they do not cross the critical line, 
solutions which start out in this region of the phase plane 
are thus candidates for spatially fiat bouncing Universes. 

Solutions below the critical line have ip < and will 
hence not exhibit a bounce. Note that the critical line 
is not itself a trajectory of the dynamics. Indeed, on the 
critical line trajectories point in vertical direction, since 
dip/dH = 0. 

There is a separatrix line between phase plane trajec- 
tories which start near the origin and which reach = 2 
(and arc thus candidates for a bouncing Universe) and 
those trajectories which cross the critical line ipc{H) and 
turn around, i.e. become solutions with -0 < 0. To 
determine the location of the separatrix line, we solve 
the equations (^ near the origin of the phase plane for 
\H\ <^ ijj, in which case the variational equation with 
respect to ip in (R) becomes 



HtP 



(18) 



which must be solved together with (jl^)- After differen- 
tiating (17) with respect to time and substituting (^8|) to 
eliminate ip, we obtain the equation 



3 



tq 




FIG. 1. The phase portrait for solutions of the equations 
of motion, with the dilaton kinetic energy (x) set to zero is 
shown. The Hubble parameter, H , is plotted on the vertical 
axis, while i/; is plotted on the horizontal. 



HH - -H^ = 0, 



which has the solution 



H{t) = -ct^ 



(19) 



(20) 



where c is a constant which labels the different trajec- 
tories. The second solution of (19) is H{t) = const and 
is thus uninteresting. From ([17D it follows that for the 
above solution 



V'(t) = ct. 



(21) 



Equations ( |20| ) and ( pi| ) give the solutions of the dynami- 
cal equations near the origin of the phase plane, provided 
that <C ■0- To get an idea for where the separatrix 
line lies, we demand that the trajectory be above the 
critical line at t/' = 2. This condition becomes 



\H\ = -?/'^ withe > 4V2. 



(22) 



Setting c = 4a/2 in ( p^ gives a first estimate for the lo- 
cation of the separatrix line. Since the set of phase plane 
trajectories near the origin is labeled by the parameter c, 
it follows from (^2|) that the set of contracting Universes 
starting near Minkowski space-time which are candidates 
for a spatially flat bouncing Universe has finite measure. 

When if) = i2, H changes sign. Contracting solutions 
which begin near the origin in the phase plane reach a 
maximal value of at "0 = 2. What happens to these 



trajectories next depends on the value of \H\ for ip = 2. 
If 

|i?(0 = 2)|«l, (23) 

then ([l^ ) gives a good approximation to the dynamics in 
the region 2 < ip < A, and we conclude that the trajec- 
tories bounce. Making use of \H\ ~ ^'^ /c (see (|2^)), the 
condition (^3|) becomes 



(24) 



O 4. 



which is consistent with the previous condition ( |22[ ) for 
bouncing trajectories. 

Trajectories with \H{'4> — 2)\ ^ 1 also reach their max- 
imal value of a.t ip — 2, but they do not bounce be- 
cause the "0 equation of motion can now be approximated 

by 



(25) 



which implies that ip keeps growing indefinitely. In com- 
bination with the equation of motion for which yields 
that 



H ~ for 101 > 4, 



(26) 



we see that the solutions tend to contracting de Sitter 
ones. Solutions below the critical line have ip < and 
will not lead to a bouncing Universe. 

For values of ip approaching ip — 2, the two critical 
lines diverge to \H\ — s- oo. This can be seen immediately 
from the ip equation of motion in (^) . For large values of 
\H\, the right hand side of this equation is dominated by 
the first two terms. In the absence of the third term, the 
solution of = which determines the critical line would 
be = 2. The third term, however, leads to a small but 
positive correction to the first term (in the bottom right 
quadrant of the phase plane which we are considering 
throughout), thus shifting the critical line slightly to the 
left of = 2, satisfying 



3H^P + 6H - —V 
H 



0. 



(27) 



To prove that all solutions of the equations (^ are non- 
singular, we must show that the solutions asymptotically 
tend to either the origin of the phase plane (Minkowski 
space-time), or to de Sitter space (the regions |0| — > oo). 
Solutions which have — > cx) for finite |0| must be 
excluded. The regions along the critical lines are the 
only part of the phase plane where there is the danger 
of singular trajectories. However, differentiating the tp 
equation of motion with respect to time, evaluating the 
result along the critical line and making use of (|2^), we 
find 



0; = -V'^ < 0. 



(28) 
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Hence, the critical line is not an attractor, but rather tra- 
jectories peel away from the line and tend to the asymp- 
totic de Sitter region. 

We have also solved the equations of motion (Q) numer- 
ically. The resulting phase diagram is shown in Figure 
1. The absence of singular solutions is manifest. In the 
asymptotic regions |?/'| ^ 4, all solutions tend to de Sit- 
ter space. The critical lines are seen to repel contracting 
solutions towards the asymptotic de Sitter regions. The 
most interesting class of solutions are those correspond- 
ing to a cosmological bounce. As predicted, they form 
a set of finite measure among solutions which start out 
close to the origin of the phase plane. There are also 
solutions which "oscillate" about Minkowski space-time. 
These solutions are further discussed in the following sec- 
tion, since we expect that they will be strongly perturbed 
by the presence of matter, in particular by the dilaton. 



IV. EFFECTS OF THE DILATON ON THE 
PHASE SPACE TRAJECTORIES 

In the presence of the dilaton, the phase space becomes 
three dimensional {ip{t), H{t), x{t)) and therefore more 
difficult to discuss analytically. We first note that in 
the Einstein frame, the dilaton corresponds to a homo- 
geneous free massless scalar field with equation of state 



P 



(29) 



where p and p denote pressure and energy density, re- 
spectively. This implies 

pit) cx a{t)-^ (30) 

or A'' = 6 in the notation of |3^. Thus, it is already 
clear from the analysis of that the dilaton will not 
introduce any singularities into the system. In fact, from 
the equation of motion for H (see (^) it follows that, 
in the region of large \H\, the presence of x will not 
change the phase space trajectories projected onto the 
(Tjj_,H) plane (a property called "asymptotic freedom" in 
pO|). However, the presence of x will greatly accelerate 
the time evolution of -0 on the given {ip,H) trajectory. 
This is easy to see for the contracting de Sitter solutions, 
since in this case the x equation of motion (see (|^) leads 
to exponential growth of x which, inserted into the tp 
equation of motion, demonstrates that at large values of 
li/^l, the x^ term dominates the evolution of ■(/'• 

We next study the effect of x on the trajectories in 
the (■0, H) phase plane. The role of the lines ip = ±2 
remains unchanged: they correspond to maxima of \H\ 
for any given trajectory. However, the condition for the 
bounce changes. Instead of V{ipb) = 0, it now follows 
from (|^) that the condition becomes 




MzL t 






+ ViiJb) = 0. 



(31) 



FIG. 2. A specific bouncing solution with a non-trivial con- 
tribution from the dilaton is plotted. Panel a) shows the evo- 
lution of // as a function of time. Panel b) depicts that of 
a(t). Panel c) that of and Panel d) that oi'^it). 
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Hence, \ipb\ is shifted to a larger (and dilaton-dependent) 
value. Note that for very large values of |x| for which the 
above equation has no solutions there will be no bounce. 

The presence of x a-lso changes the critical lines. For 
a fixed value of Tp, the condition ^ — which determines 
the critical line occurs at a larger value of than it 
does in the absence of x, as can easily be seen from the iJj 
equation of motion (see (^. In addition, for fixed initial 
{ip,H) the value of ^ for contracting solutions is larger 
with X 7^ than with x = 0, so we conclude that adding 
a small value of x increases the range of initial conditions 
in the H) plane which lead to a bounce. Furthermore 
the larger the value of x, the larger the effect is (as long 
as ( |3l| ) still has a solution). 

We will now argue that collapsing solutions in the pres- 
ence of a small |x| quite generically lead to a spatially 
flat bouncing Universe. Consider initial conditions with 
small but positive x which in the H) plane lie below 
the separatrix line and which initially "oscillate" about 
Minkowski space-time. Note that since 



(32) 



and since -0 > for these trajectories, they do not bounce 
upon reaching H — but start another cycle with H < 0. 
Since < at all times, x{t) is increasing. Eventually, 
x(t) reaches a sufficiently large value such that the tra- 
jectory crosses the "separatrix" in the ("0, H) plane and 
evolves past = 2 to a successful bounce. An example 
of such a trajectory is shown in Figure 2. 

It is now obviously possible to construct initial con- 
ditions in which the Universe is initially dominated by 
the dilaton and contracting towards a bounce. For ex- 
ample, we can take the initial conditions for (■0, TJ, x) to 
be those of the solution in Figure 2 at the end of the cy- 
cle preceding the bounce. We have thus shown that our 
model provides a successful implementation of the evo- 
lution postulated in pre-big-bang cosmology. Note that 
after the bounce, the dilaton rapidly becomes irrelevant 
to the evolution in the expanding phase. Note, in partic- 
ular, that the dilaton tends to a constant. This feature is 
different from what happens in other attempts to achieve 
a branch change where the dilaton continues to grow in 
the expanding phase. 



V. CONCLUSIONS 

We have shown that the "limiting curvature construc- 
tion" of Rcfs. 30 1 can be applied to dilaton cosmology 
to yield spatially flat bouncing cosmological solutions. In 
particular, this enables us to implement a model which 
successfully implement the qualitative evolution associ- 
ated with a branch change in pre-big-bang cosmology. 

Speciflcally, we have investigated spatially homoge- 
neous, isotropic, and flat solutions of the equations of 
motion for dilaton gravity coupled to, I2, a special com- 
bination of invariants quadratic in the Riemann curva- 



ture. The coupling is provided by a Lagrange multiplier 
field ip. We have determined the conditions which must 
be imposed on the potential V{il)) required in order to 
guarantee that all solutions are singularity-free, and to 
permit bouncing scenarios. 

We then studied the phase space of trajectories and 
demonstrated, both analytically and numerically, that 
all the solutions to the equations of motion which follow 
from our action are nonsingular. We identified a large 
class of contracting trajectories which start out near the 
origin in phase space and develop a cosmological bounce. 
The addition of the dilaton increases the fraction of phase 
space for which solutions bounce. These bouncing cos- 
mologies display the qualitative behavior associated with 
the pre-big bang branch change. 

An interesting result of our work is the conclusion 
that the dilaton is not necessary to obtain a spatially 
flat bouncing cosmology. Also, in many cases a dilaton- 
dominated bounce is preceded by a period when the dila- 
ton has a negligible effect on the dynamics of tp{t) and 
H{t). 

Our work does not address the issue flxing the value 
of the dilaton at late times. Nor does it address recent 
objections |pl|-|33[| to the pre-big-bang scenario. In work 
in progress, we are investigating the possibility that our 
implementation of the branch change will alleviate some 
of these difficulties, in particular the flatness problem 
and the problem of how to obtain a large Universe. The 
relevant feature of our model which we expect to play an 
important role in addressing these problems is that the 
dynamics of our bounce are determined by the higher 
derivative gravity terms which lead to de Sitter phases 
for many trajectories. Thus, on balance, our model looks 
quite promising. 

The obvious drawback of our method is that the extra 
terms in the action are put in by hand rather than being 
derived from fundamental physics. There are, however, 
several justiflcations for following our approach. First, 
it is generally expected that physical invariants must be 
limited in string theory. This should, at high curvatures, 
be reflected in the effective action for gravity, and our 
action is an easy way to realize limiting quantities. Sec- 
ondly, higher derivative terms will inevitably arise in ef- 
fective theories of gravity, such as those derived from 
string theory, from other approaches to quantum grav- 
ity, or from quantizing matter fields in curved space-time. 
Our action, even though not directly derived from string 
theory, demonstrates that a specific set of higher curva- 
ture corrections can both ensure that all relevant physical 
invariants remain finite, and produce a bouncing nonsin- 
gular cosmology. 
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